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SUMMARY 

Development of probabilistic structural analysis methods for hot structures is 
a major activity at NASA Lewis Research Center. It consists of five program ele- 
ments: (1) probabilistic loads, (2) probabilistic finite element analysis, and (3) 

probabilistic material behavior, (4) assessment of reliability and risk, and (5) 
probabilistic structural performance evaluation. Recent progress includes: (1) 

quantification of the effects of uncertainties for several variables on High Pressure 
Fuel Turbopump (HPFT) blade temperature, pressure, and torque of the Space Shuttle 
Main Engine (SSME), (2) the evaluation of the cumulative distribution function for 
various structural response variables based on assumed uncertainties in primitive 
structural variables, (3) evaluation of the failure probability, (4) reliability and 
risk-cost assessment, and (5) an outline of an emerging approach for eventual hot 
structures certification. Collectively, the results demonstrate that the structural 
durability/reliability of hot structural components can be effectively evaluated in a 
formal probabilistic framework. In addition, the approach can be readily extended to 
computationally simulate certification of hot structures for aerospace environments. 


INTRODUCTION 

It is becoming increasingly evident that deterministic structural analysis 
methods will not be sufficient to properly design critical components in hot struc- 
tures in general and propulsion structures in particular. These structural com- 
ponents are subjected to a variety of complex, and severe cyclic loading conditions, 
including high temperatures and high temperature gradients. Most of these are 
quantifiable only as best engineering estimates. These complex loading conditions 
subject the material to coupled nonlinear behavior which depends on stress, temper- 
ature, and time. Coupled nonlinear material behavior is nonuniform, is very dif- 
ficult to determine experimentally, and it is difficult (if not impossible) to 
quantify deterministically. In addition, hot rotating structural components are 
relatively small. Fabrication tolerances on these components, which in essence are 
small thickness variations, can have significant effects on the component structural 
response. Fabrication tolerances by their very nature are statistical. Furthermore, 
the attachment of components in the structural system generally differs by some 
indeterminant degree from that assumed for designing the component. In summary, all 
four fundamental aspects — (1) loading conditions, (2) materials behavior, (3) 
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geometric configuration, and (4) boundary conditions (on which structural analyses 
are based) — are of random nature. The direct way to formally account for all 
these uncertain aspects is to develop probabilistic structural analysis methods where 
all participating variables are described by appropriate probability functions. 

The development of the probabilistic structural analysis methodology is an 
on-going joint program of NASA Lewis Research Center in-house and sponsored research 
(ref. 1). Theoretical considerations, computer codes, and other relevant applica- 
tions are described in papers presented in various recent conferences (refs. 2 to 
32). Activities and progress up to June 1989 are summarized in reference 33. The 
objectives of this invited paper are (1) to provide a brief description of the fund- 
amental aspects and (2) to demonstrate the application of this methodology to a 
specific example (the reliability/risk of turbine blade components of rocket pro- 
pulsion systems). The specific example includes the four key elements required in 
probabilistic structural analysis of thermal structures, namely: (1) probabilistic 

loads, (2) probabilistic finite element analysis, (3) probabilistic description of 
coupled nonlinear material behavior, and (4) evaluation of reliability and risk. 
Throughout the specific example discussion appropriate comments are included to 
illustrate the generality of the method and its application to a variety of hot 
structures . 


FUNDAMENTAL CONSIDERATION 

Central to the probabilistic structural analysis described herein is the funda- 
mental consideration that: Uncertainties observed in the structural performance 

(displacements), frequencies, buckling, global fracture toughness, stresses/strains 
of hot structures can be quantified in terms of corresponding, uncertainties in basic 
parameters (primitive variables) . The primitive variables are those which are used 
to describe the structure and its respective environment. For example: (1) struc- 

tural configuration, (2) boundary conditions, (3) loading conditions, and (4) 
material thermomechanical nonlinear behavior. 

The uncertainties in these primitive variables are then integrated through 
structural mechanics to quantify the uncertainties in the global structural responses 
(example-displacements) and are decomposed to quantify the uncertainties in local 
responses (example-stresses/strains) . The concept is schematically illustrated in 
figure 1. The structural component is the blade which is modelled for finite element 
analysis. The input uncertainties are the blade loads (centrifugal, pressure, and 
temperature), geometry and material variables. The output is quantification of 
uncertainties in structural responses or in local stresses for probable fracture 
initiation. Brief descriptions are given on subsequent sections on each of these as 
it is applied to the specific example. 


PROBABILISTIC SIMULATION OF LOADS 

The fundamental assumption for the probabilistic simulation of loads is that 
each individual load condition is the probabilistic time synthesis of four primitive 
parts: (1) steady state, (2) periodic, (3) random, and (4) spike. Each of these 

parts, except random, is described by a deterministic portion and a probabilistic 
perturbation about this deterministic portion. The resulting distribution is similar 
to the schematic, in figure 1 upper left and as described in reference 33. One 
justification for describing each loading condition in terms of primitive parts is 
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that experts, over the years, have developed good judgment on the ranges of the 
perturbations about nominal (deterministic) conditions. A computer code with 
dedicated expert systems (Composite Load Spectra) has been developed to synthesize 
these four parts by using (1) available data from various SSME engines, (2) probabil- 
ity theory, and (3) expert opinion, as depicted schematically in figure 2. 

The results from the application of the Composite Load Spectra computer code to 
probabilistically simulate loads for two blades are summarized in table I. The first 
5 lines are uncertainties in some of the engine specific primitive variables, while 
the last four are the predicted loads. The comparisons of predicted loads with the 
measured data are in very good agreement considering the large number of primitive 
variables (47) required to synthesize these loads. The conclusion is that methods 
can be developed and are available to probabilistically synthesize complex load 
conditions for hot structures. 


PROBABILISTIC FINITE ELEMENT STRUCTURAL ANALYSIS 

The fundamental assumption for developing probabilistic finite element 
structural analysis ( PFEM) is that the uncertainties in each primitive structural 
variable can be represented by a probabilistic distribution. Primitive structural 
variables are those which are used to describe a structure such as: (1) stiffness, 

(2) strength, (3) thickness and tolerance, (4) spatial location, (5) attachment, and 
(6) various nonlinear material dependencies (temperature, stress, time, etc) typical 
to the schematic in figure 1 upper right (ref. 33). Subsequently, the uncertainties 
in the load conditions (described by the composite load spectra) and the uncer- 
tainties in the primitive structural variables are computationally synthesized by 
performing probabilistic finite element structural analysis to simulate uncertainties 
in the structural response of a specific SSME structural component. The structural 
response such as displacement, frequencies, buckling loads, and global fracture 
toughness as was already mentioned, is generally described in terms of cumulative 
probability distribution. 

PFEM has been formalized and integrated into a computer code identified as 
NESSUS (Numerical Evaluation of Stochastic Structures Under Stress). NESSUS is 
driven by an expert system. A schematic diagram of NESSUS is shown in figure 3. 
Representative results obtained using NESSUS for one of the blades are shown in 
figure 4 for stress at two different points. The distributions assumed for the 
primitive variables are listed in table II in addition to the blade loads from the 
CLS computer code. The combined effects of the primitive variable uncertainties on 
effective effect stress are shown in figure 4 in terms of cumulative probability 
distribution (CPD). The information generated during the generation of the CPD can 
also be used to evaluate the sensitivities as shown in figure 4 for each point. The 
significant point is that the PFEM yields a wealth of information which can be used 
to evaluate: (1) the uncertainties in the local structural response, and (2) the 

sensitivities which can be used to adjust the design for enhanced probability of 
success. The important conclusion is that probabilistic finite element methods can 
be developed and are available to quantify uncertainties in the structural perfor- 
mance of hot structures. In addition the sensitivities in the primitive variables 
that influence this performance can be evaluated and ranked. 
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PROBABILISTIC SIMULATION OF MATERIAL BEHAVIOR 


The fundamental assumptions to probabilistically simulate nonlinear material 
behavior ares (1) a relationship for material behavior can be developed in terms of 
primitive variables affecting this behavior and (2) the uncertainties in the primi- 
tive variables can be described by assumed distributions (ref. 33). A multifactor 
interaction model (MFIM) for this relationship is shown in figure 5, where its con- 
stituent primitive variables are identified. 

This MFIM is applied to the blade to probabilistically determine: (1) the 

resistance curve for damage (crack) initiation and (2) damage propagation and its 
effects on global structural response. The inputs are summarized in table III. The 
results for the most probable damage path, with a probability of occurrence of 

0.0002, and the respective degradation in frequencies are shown in figure 6. The 
important observation is that the uncertainties in damage initiation, propagation and 
subsequent effects on structural performance can be probabilistically simulated by 
the methodology described herein. It is worthy of note that this methodology can 
readily be incorporated to monitor the in-service health of aerospace hot structures 
(provide suitable monitoring devices are available). 


RISK-COST ASSESSMENT 

The methodology described previously can be extended to perform reliability and 
risk-cost assessments." in order to accomplish this, (1) the cost for component 
flight readiness needs to be quantified and (2) the cost as a consequence of failure 
must be established. Assumed aspects of both of these have been integrated into the 
methodology described herein (ref. 34). Application of this methodology to the blade 
in the previous section yields the results summarized in figure 7. 

It can be seen in figure 7 that: (1) the probability of damage initiation can 

be evaluated versus fatigue cycles and (2) the total cost, which is used to assess 
risk, is evaluated versus fatigue cycles. The important observation from the afore- 
mentioned discussion is that the reliability and risk-cost of hot structures can be 
assessed using the probabilistic methods of the type described herein. The implic- 
ations are far-reaching because these methods can be applied to existing hot struc- 
tures to evaluate their risk for continuing service as well as those on the design 
board and those still in the conceptual phase. 

RELIABILITY/CERTIFICATION — AN EMERGING APPROACH 

The collective observations from the previous discussion lead to an emerging 
approach to computationally simulate hot structures reliability, risk, component 
qualification, and eventually vehicle structure certification. The general steps for 
this emerging approach are outlined as follows: 

1. Develop a hot structural component/vehicle (global) analysis model. 

2. Conduct probabilistic structural analysis (PSA) of the types described 
herein. 

3. Identify the critical component/vehicle areas from the results of PSA. 
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4. Perform global/local PSA to evaluate nonlinear effects and to locate 
probable sites of damage initiation, 

5. Determine the most probable damage propagation path, 

6. Evaluate probable structural degradation along this path, 

7. Establish probable path extent for violation of specified structural per- 
formance criteria (for example, 10-percent reduction in the first vibration mode), 

8. Assess corresponding reliability and risk and decide on their acceptability. 

9. Schedule inspection intervals and retirement for cause criteria based on the 
results of items 5 to 7 , 

10. Verify with probabilistically selected (using respective sensitivities) 
critical hot-structures components and prototype tests. 

11. Design a suitable in-service health monitoring system using the results 
from items 8 and 9 above in order to ascertain that the component/structure will meet 
the acceptable reliability and risk. 


CONCLUSIONS 

A methodology has been developed for the formal probabilistic quantification 
of uncertainties in the structural performance of aerospace hot structures. The key 
elements in this methodology are: (1) probabilistic load simulation, (2) probabil- 

istic finite element analysis, (3) probabilistic simulation of thermomechanical 
nonlinear material behavior, and (4) risk-cost assessment. This methodology is 
described herein in terms of fundamental aspects and application to a specific struc- 
tural component which is a turbopump blade of the Space Shuttle Main Engine (SSME) . 
The specific example illustrates how the uncertainties in all the basic parameters 
(primitive variables) for loads, structure and material behavior are incorporated in 
order to probabilistically simulate the uncertainties in the structural response 
(global and local). Also, the example illustrates how the reliability and risk-cost 
can be assessed. Collectively, the summary of the fundamental considerations and the 
results from the specific example demonstrate that a formal methodology is available 
to evaluate the reliability and risk-cost of hot structures in aerospace environ- 
ments. In addition, an emerging approach is outlined which can be used to computa- 
tionally qualify and eventually certify hot structures. 
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TABLE I. - HPOTP AND HPFTP PARAMETERS PHASE II ENGINE CALCULATED VERSUS MEASURED 


Condition 

HPOTP 

HPFTP 


Speed 

Turbine 

Speed 

Turbine dis- 




discharge 



charge temper- 




temperature 

- 


ature 


Calcu- 

Meas- 

Calcu- 

Meas- 

Calcu- 

Meas- 

Calcu- 

Meas- 


lated, 

ured. 

lated, 

ured, 

lated. 

ured, 

lated, 

ured, 


rpm 

rpm 

°R 

°R 

rpm 

rpm 

°R 

°R 

Hardware - 20 random 

294 


53 

— 

288 


65 

— 

Test - 20 random 

210 


157 


306 


20 

— 

Total random 

360 


165 

— 

554 


70 

— 

Low NPSP - det 

620 


225 

— 

56 


52 

— 

High NPSP - det 

-317 


-219 


-94 


-62 

— 

Range « random + det 

1 660 

1 500 

475 

400 

1 260 

1 000 

114 

150 

Max 

28 090 

29 100 

1630 

1650 

35 742 

35 750 

1740 

1760 

Norn 

28 100 

28 200 

1374 

1380 

35 130 

35 300 

1688 

1690 

Min 

27 430 

27 500 

1155 

1250 

34 482 

34 750 

1625 

1610 


Measured: Measured variation for phase II test set. 

Hardware: Variations engine hardware. 

Test: Initial lost conditions - inlet temperatures and mixture ratio. 

Det: Duty cycle effects of inlet pressures plus correlated 2 O variations of 

cavitation. 
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TABLE II. - RANDOM VARIABLES CONSIDERED AND THEIR STATISTICS 


Number 

Random variable 

Type 

Affected FEM quan- 
tities 

Mean 

Standard 
deviat ion 

1 

Material axis Z 

Material 

Anisotropic 

-0.087266 radian 

0.067544 

2 

Material axis Y 

orientation 

material 

-0.034907 

0.067544 

3 

Material axis X 

Effects 

Orientation angles 

-0.052360 

0.067544 

4 

Elastic modulus 

Material 

Elastic constants 

18.38E6 psi 

0.4595E6 

5 

Poisson's ratio 

properties 


0. 386 

0.00965 

6 

Shear modulus 



18 . 63E6 psi 

0 . 046575E6 

7 

Geometric lean 

Geometrical 

Node coordinates 

0 deg 

0.14 deg 

8 

Geometric tilt 

variations 


0 deg 

0.14 deg 

9 

Geometric twist 



0 deg 

0.30 deg 

10 

Mixture ratio 

System 

Pressure, temperature. 

6.0 

0.02 

11 

Fuel inlet pressure 

independent 

centrifugal force 

30.0 psi 

5.00 

12 

Oxidizer inlet pressure 

loads 


100.00 psi 

! 26.00 

13 

Fuel inlet temperture 



37 ° R 

! 0.50 

14 

Oxidizer inlet temperature 



164 ° R 

1.33 

15 

Pump efficiency 

' Component 

Pressure, temperature, 

1.00 

0.008 

16 

Head coefficient 

independent 

centrifugal force 

1.024 

0.008 



loads 




17 

Coolant seal leakage 

Local effects 

Temperature 

1.0 

0.10 

18 

Hot gas seal leakage 



1.0 

0.05 


TABLE III. - PRIMITIVE VARIABLE PROBABILITY DISTRIBUTIONS 
FOR PROBABILISTIC MATERIAL PROPERTY MODEL 


Variable 

Distribution 

type 

Mean 

Standard 

deviation 

Value 

Percent 
of mean 

T 

F 

Normal 

2750 °F 

b 

o 

<4* 

in 

2.0 

T 

o 

Normal 

CL 

o 

GO 

VD 

fO 

o 

o 

3.0 

S 

p 

Normal 

212.0 ksi 

10.6 ksi 

5.0 


Constant 

0 

0 

0 

N 

MF 

Lognormal 

10 8 

5xl0 6 

5.0 

N 

MO 

Lognormal 

10 3 

50 

5.0 

n 

Normal 

0.25 


3.0 

P 

Normal 

0.25 


3.0 

q 

Normal 

0.25 


3.0 
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Turbine blade loading 


Nessus turbine blade coarse model 


Geometry and 
material behavior model 





Resistance strength - S 
(from generic 
probabilistic 
aterial behavior model) 


Response s 

L Probable 

loan 


amage initiation! 


LeRC contacts: 

CLS - Composite Load Spectra - Rocketdyne 

PSAM - Probabilistic Structural Analysts Methods - SWRI 

PMBM - Initial assessment - In-house 


|-« — Resistance scatter— ►] 


Figure 1. — Probabilistic simulation of component reliability using CLS coupled with PSAM and PMBM. 


Composite load spectra expert system 


Expert 

system 



LDEXPT design 

Load expert system design philosophy 

• Rule-base production system 

- If then rules 

• Simple inference scheme 

- Inference net (decision tree) 

• Sophisticated probabilistic methods 

- Discrete probability distribution 

- Monte Carlo 

- Barrier crossing 

« Powerful knowledge base 

- influence coefficient 

- Scaling coefficients 

- Duty cycle load profiles 

- Engine configuration and geometry data 

- Raw engine flight and test data 


Figure 2. — Composite load spectra simulation using expert systems. 
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Cumulative probability 



Figure 3.— Computational simulation of probabilistic 
structural response using Nessus. 


Turbine blade loading Nessus turbine blade coarse model 


Geometry and 



Material axis 0 y 0.624 

Geometric lean angle .478 
Material Axis 8 Z .320 
Fuel inlet temperature .220 



Hot gas seal leakage 0. 776 
Material axis 0 y .349 

Material axis 6 Z .298 

Poisson's ratio .247 



Effective stress, ksi 


Effective stress, ksi 


Figure 4. —Probabilistic component stress analysis. 
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Probability of occurrence Probability of occurrence 


Mp = Mpo 


1 

-n* 

1 

H 

n 

Sf — G 

p " 

LTf-ToJ 


o 

[to 



logN MF -k>gN M 
log Hjif -k>g Nmq 


Primitive variables Subscripts 

Mp = Degraded material property F = Final characteristic value 

T = Temperature O = Reference property 

S = Strength 
a = Stress 
N m = Mechanical cycles 

Figure 5. — Generic probabilistic material property degradation 
model in terms of primitive variables. 





Natural frequencies, Hz Natural frequencies, Hz 


Figure 6. — Natural frequencies decrease as fracture progresses. 




Damaged 


state 


(undamaged) 

O 

D° 

□ 

D 10 

A 

D 10, 9 

V 

Qto, 9, 14 

O 

plO, 9, 14, 18 
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Probability of 
damage initiation 

o o, o 

U ro 

— L L 1 

10 4 10 5 10 6 
Fatigue cycles 

Total cost 

= Cost for components 
service readiness 

+ Probability of damage 
initiation 

* Consequential cost 
due to damage 
initiation 





Fatigue cycles 


Figure 7. — Probabilistic risk-cost assessment. 
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